Localization and derealization in the quantum kicked prime number rotator 
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The quantum kicked prime number rotator (QKPR) is denned as the rotator whose energy levels 
are prime numbers. The long time behavior is decided by the kick period r and kick strength k. 
When is irrational, QKPR is localized because of the equidistribution theorem. When ^- is 
rational, QKPR is localized for small k, because the system seems like a generalized kicked dimer 
model. We argue for rational ^ QKPR delocalizes for large k. 
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The kicked prime number rotator is denned as 



H = H -V^T6{t-nT), 



(1) 



where .Ho is the unperturbed Hamiltonian, and V is 
the perturbation. Hq is a diagonal matrix. The ru- 
th eigenvalue E m corresponding to m-th eigenstate \m) 
of Hq is the |m|-th prime number p\ m \. When m < 
0, E m = E- rn . E = 0. The diagonal of H is 
{...,11,7,5,3,2,0,2,3,5,7,11,...}. 
V is defined as 
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The Floquet operator is F — e n v ( e ) e n H ° T . The ma- 
trix elements of F is F nm = exp(=^)i m -™ J n _ m (£), 
where E m — p\ m \, J n -m is the Bessel function of the first 
kind. We set % = 1. F nrn — exp(— iTE m )i m ~ n J„_ m (fc). 
The system is very like the quantum kicked rotator 
(QKR), except its energy levels are now prime numbers. 

It seems there is no classical correspondence of QKPR. 
Experimental implementation of such a model also seems 
impossible. Nevertheless it still has some theoretical in- 
terests. In the paper, we numerically calculate the evo- 
lution of QKPR. We are interested in the same problem 
in QKR. If the particle is in the ground state |0) initially, 
will it diffuse away in the future? 

The evolution of the system is calculated by the iter- 
ative unitary matrix multiply method [l]. F 4 = (F 2 ) 2 . 

F 4 = ^2)2_ ^8 = ( F 4)2^ And SQ on ' In thig way; we 

can calculate F 1 - 2 ^ ' by 50 matrix multiplies. In all our 
calculation, N indicates at time 2 n t. For example, the 
first figure N = 10 means at time 2 10 t = 1024r. n is 
the n-th basis \n) and c„ is the base-10 logarithm of the 
absolute value of the wave function on the \n). n runs 
from —500 to 500 in our calculation. 
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The result is 



First, we choose k = 1, r = 2tt 
displayed in FIG. 1. QKPR is localized perfectly. In our 
simulation, the exponentially fall of the wave function 
never changes from N = 8 to N = 50. The wave function 
on the |n) is 10 c ™. From iV = 1 to N = 7, the wave 
function is somewhat curved. After the first kick, the 
wave function is the 0-th column of the Floquet matrix 
F. The absolute value of F n0 is |J n (fc)|. \J n (k)\ falls to 
zero faster than exponentially . This is the reason the 
curved form of the wave function. 

Second, we choose k = 5, r — 2n ^' 2 ~ 1 . The result is 
displayed in FIG. 3. The wave function is also localized. 
This is expected. The sequence {— p n j^Modl} is cquidis- 
tributed between [0, 1], when ^ is irrational. We denote 

the sequence {-pn^j^-Modl} as QKPR G . In QKR, 

the sequence {-j^Modl} is also equidistributed be- 
tween [0, 1], for an irrational We denote the sequence 

{-^-^^Modl} as (QKR G ). We can also use the in- 
verse Cayley transform method to convert the Floquet 
eigenstate equation Ftp = X(p into an equation like An- 
derson localization problem. From Fishman et aVs argu- 
ment 0, QKPR will localize. 

In the left of FIG. 2, QKPR G and QKR G are dis- 
played. Though there are apparently some correlations 
in QKPRq and QKR G and the correlation is different 
between both sequences. The correlation is surely not 
strong enough to destroy localization. If a sequence is 
periodic with a period q, then the discrete Fourier trans- 
form of the sequence is composed by q modes. To find 
whether there is some periodicity in the sequence, we 
perform a discrete Fourier transform (DFT) on the se- 
quence. DFT of a sequence s n of length L is defined as 



-i27r(n-l)(j-l)/L 



where j runs from 1 to 
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L. There are some other definitions of DFT with nu- 
anced difference with our definition. But the difference 
is irrelevant to our discussion here. In the right of FIG. 
2, the DFTs of both sequences are displayed. There are 
no rigorous periodicity in both sequences. F^ of QKPRq 
seems to have a trend to cluster together. Also it is less 
uniformly distributed than the Fj of QKRg and tends to 
be small. 

If = | , does QKPR localize? At first thought, 
this seems to be a resonant case in QKR and the rotator 
will delocalize. The calculation result is in fact it still 
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localizes for small fc. In FIG. 4, we choose k = 1 and 
r = 2tt± and in FIG. 5, k = 5 and r = 2tt±. QKPR of 
fc = 1 is apparently localized. 

The explanation of the localization is QKPR with 
= 3 and fc = 1 is a kicked pseudo dimer rotator. 
The dimer model is defined as every diagonal matrix ele- 
ment is a probability variable which only takes two values 
[1, H[ . The kicked dimer model can be defined as every 
diagonal matrix element of e~ tH ° T is a random variable 
which takes two values. If it takes more than two values, 
it is a generalized kicked dimer model. For q — 5, the 
sequence — p n ^Modl mainly takes four values. So it is 
a generalized kicked pseudo dimer model. The sequence 
—p n TModl (QKPR3) is not really random. But thepseu- 
dorandomness is enough to result in localization [2j]. To 
measure how random QKPR3 is, we perform a DFT on it. 
In the left of FIG. 4, we compare QKPR3 with a dimer 
sequence D3, which is defined as —^Random(n)Modl, 
where every Random(n) is a random variable which takes 
two values 1 and 2. The F k s of D3 and QKPR3 are quite 
close with each other, except QKPR3 tends to cluster 
together. 

Does QKPR localize for ^ = | and fc = 5? We 
think it delocalizes. There are a series of plateaux in 
the wave function of QKPR. The wave function falls 
abruptly when approaching the boundary (cliff) of a 
plateau. Some plateaux disappear at N — 14,27,45. 
QKPR wave pass through the cliff, so it disappears in- 
termittently. 

The most obvious cliff is from n = 140 to 160 in FIG. 6. 
QKPR 3 from n = 139 to 161 is 
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Note it is periodic from n — 140 to 160. In a perfect 



periodic potential the wave will always diffuse away. The 
quantum wave can not stay at a period potential very 
long. Once the wave has propagated into the phase space 
between n = 140 and 160, it diffuses away quickly. While 
once the wave propagate the phase space whose neigh- 
borhood —pn/'&Modl is irregular, the quantum wave is 
localized there at least temporarily. When k = 1, n from 
140 to 160 is a plateau in FIG. 2. There is a transition 
from a plateau to a cliff when fc becomes large. 

Note in FIG. 4, for lots of fc, Lots of F k of QKPR 3 and 
QKPRg tend to be small. From Plancherel's theorem, 
there must be some F k tends to very large. So there 
is some weak periodicity in the sequence QKPR3. The 
sequence from n = 140 to 160 increases the periodicity 
of QKPR3. 

At N = 54, there is a plateau from —50 to 50 or so. 
Even at N — 50, there is a plateau from —50 to 50. The 
localization length of QKR is ^ = 6.25 or so [1,11,0]. 
So the distribution length of QKPR is much larger than 
the localization length. Even the result of N = 54 is 
untrustworthy, we think the quantum wave is absolutely 
not localized in the localization length. From n = —150 
to 150 or so, there is apparently unneglectable quantum 
wave at every n at for example N = 40. At N = 6, 
the quantum wave has already propagated into n > 6|. 
QKPR with f = § and fc = 5 is not localized. In [g], 
we point out when the kick strength is larger than n, the 
inverse Cayley transform method breaks down. QKPR 
is apparently an evidence to the failure of the inverse 
Cayley transform method when fc is large. 

In this paper, we apply the iterative unitary matrix 
multiply method to quantum kicked prime number rota- 
tor. If is irrational, the rotator localizes. If is ratio- 
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nal, for small kick strength fc, the rotator localizes. As fc 
increases, we argue there is a localization-delocalization 
transition. 
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FIG. 1: QKPR wave function at different time for k = 1, r = 2ir^^-. 
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FIG. 2: Left: —E n r modulo 1, where E„ is the n-th energy level. Right: DFT of the sequence {exp(— iE n r)}, where n 
runs from 1 to 500. 




FIG. 5: QKPR wave function at different time for k = 5, r = 27r|. 
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FIG. 6: QKPR wave function plateaux and cliffs at N = 20 and N = 40 for k = 5, r = 27ri. 



